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We propose an implementation sclieme for tlie continuous-time quantum walk using a diatomic 
molecule and an optical frequency comb. We show an analogy between the quantum walk and 
the cascade rotational transitions induced by the optical frequency comb whose frequency peaks 
are tuned to the pure rotational transitions in the molecule. The strategy to compensate for the 
centrifugal distortion of the real molecule is also demonstrated. 
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I. INTRODUCTION 

The quantum walk is the quantum counterpart of the well-known classical random walk. The quantum walk has 
been actively studied for the quantum information processing [l|, 13]. Historically, two types of quantum walks have 
, been studied independently, that is, the continuous-time quantum walk (CTQW)^ and the discrete-time quantum 
Q-c walk (DTQW) [j] . The CTQW describes a continuous diffusion of the probability amplitude on a series of discrete 
sites, whereas the DTQW describes an expansion of the entanglement by the discrete unitary operations on a series 
^ of qubits. In these days, the concept of the CTQW attracts more and more attention. For example, the CTQW can 
be regarded as a universal computational primitive i5]. Also, the concept of the CTQW helps us to understand deeply 
the energy transfer mechanism in photosynthesis @, |7j. 

The implementation schemes of the CTQW were studied by some groups. The 4-state cyclic CTQW was realized 
first on a two-qubit nuclear- magnetic- resonance quantum computer [8]. Schemes based on the quantum dotsj^, and 
Rydberg atoms in an optical lattices jlOil were also suggested. Recently, the large scale one-dimensional CTQW was 
_ realized by using waveguide lattices The waveguide lattices can be used for the advanced information processing 
ly-j ■ using correlated photons [l^. The implementation of the quantum walk should facilitate the experimental simulation 
^-H ' of various phenomena, even if that cannot be used for the correlated particles. To be used for the simulation system, 
. both the scalability and the controllability are required for the implementation scheme. To determine the probability 
.^-^ ' distribution of the CTQW of many steps, parallel execution of a huge number of trials is favorable because the number 
, of required trial rapidly increases as the number of steps increases. Also, to implement various modifications, the 
T-H ' driving force of the CTQW should be controllable. In this paper, we propose a new implementation scheme of the 
^-H ] CTQW using a diatomic molecule and an optical frequency comb. Both the controllability and the scalability are 
^ ■ expected to be included in this scheme. 
• I— I ' 

X: 

H ■ II. THE CONTINUOUS-TIME QUANTUM WALK 

a • 

We briefly review the definition of the CTQW. In one dimensional space, the CTQW is expressed as the finite- 
difference Schrodinger equation 



(1) 



where t) is a complex amplitude at the continuous time t and discrete site position n, and 7 is the rate constant 
of the transition between neighboring sites. With the initial condition of ijj{n,0) = dn,o, the exact solution of CTQW 
is given by 

V'(n,t) = z'"l:^„|(70, (2) 
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FIG. 1: (Color online) An example of probability distributions of the continuous-time classical random walk and the CTQW. 



quite differently from its classical counterpart (Fig.l). The main features of the CTQW, compared with the classical 
random walk, are the peaks that stands at both ends, the rapid diffusion of the probability distribution, and the sharp 
reflecting behavior at the boundary. 



Diatomic molecule is one of the most characterized multi-level system. First, we regard the rotational states of a 
diatomic molecule as the nodes of the quantum walk. Due to the selection rule of AJ — ±1, the rotational states 
compose a one-dimensional network which is connected by the pure rotational transitions. When the molecules are 
irradiated by the electromagnetic wave which can induce both AJ = ±1 transitions, the molecules are excited, de- 
excited, or remained as it is according to a certain probability. Furthermore, if the electromagnetic wave induces 
all the rotational transitions in the network, the irradiated molecules are repeatedly excited or de-excited, as if they 
were randomly walking on the network. However, the evolution dynamics of them is not the random walk but the 
quantum walk. To induce all the rotational transitions in the molecule, a broad-band laser pulse seems to be suitable. 
Actually, however, the broad-band laser pulse usually induces multi-photon processes such as ionization. 

To avoid the undesired processes, we consider the optical frequency comb which can induce only a series of the pure 
rotational transitions (Fig. 2). We analytically show the analogy between the CTQW and the cascade pure rotational 
transitions induced by the optical frequency comb in the molecules. At first, we discuss an ideal situation, in which the 
centrifugal distortion can be fully ignored. The rotational energy of diatomic molecules with a closed shell structure is 
expressed as Ej = hBJ{J+ 1), where B is the rotational constant. Because the transition frequency for J ^> J -I- 1 is 
expressed as i^j — 2B{J + 1), such transitions exhibit an evenly spaced spectrum like a comb in the frequency domain. 
Under the dipole approximation for linearly polarized electric field e(t), the equation of rotation in the interaction 
picture is written by 

ih^cj{t) = -^,/_i£(t)cxp(27rjz/j_it)cj_i(t) 



where cj{t) is the complex amplitude of the rotational state J, and /ij is the transition dipole moment for the pure 
rotational transition from J to J + 1, which is calculated by = fj,{ [(J -I- 1)^ — M^] / [(2 J -I- 1)(2 J + 3)]}^/^ using 
the permanent dipole moment /i and the magnetic quantum number M. In this paper, /x = 1.0, M = 0, and h = 1.0 
are assumed for simplicity. 

For the optical frequency comb, we define the electric field as 




III. PHYSICAL IMPLEMENTATION USING DIATOMIC MOLECULES 



Hje{t) eyip{-2TTWjt)cj+i{t), 



(4) 




(5) 



The e{t) is regarded as a train of pulses whose interval is (2_B) ^ in the time domain (Fig. 3). 
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FIG. 2: (Color online) Conceptual picture and schematic diagram of the implementation scheme of the CTQW using the 
diatomic molecules and the optical frequency comb in the ideal situation. 
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FIG. 3: (Color online) Time profiles of the optical frequency comb with various numbers of frequency component J, 



Using eqlH and putting all the oscillating terms into term A (rotating wave approximation) , eq|4] is re- written by 
«^cj(t) = -|c;_i(t) - |cj+i(t) - A. (6) 

The first two terms in eq|6] represent the contribution from the resonant frequencies (resonant contribution), whereas 
the third term A represents the contribution from the detuned frequencies (detuning contribution). The resonant 
contribution has the same form as the definition of the CTQW, eqlU by letting cj{t) = ip{n,t) and A J = n. 

The detuning contributions are expected to be averaged out to zero by integrating over the duration of {2B)~^ 
because all the exponential terms in A seems to oscillate with a period of {2B)~^. We confirmed this speculation 
numerically. Eq|4] and eqlU were solved by the fourth order Runge Kutta algorithm with the initial condition of 
cj{—0.5/2B) = Sj^ioo, Jmax = 200, and 7 — 2B. We have found that the simulated population distribution after the 
50-th pulse completely overlapped with the exact solution as shown in Fig. 4. Therefore, we conclude that the optical 
frequency comb implements the CTQW on the rotational states in a diatomic molecule under the ideal situation. 

Also, we briefiy discuss the contribution of the two-photon transition, which has been ignored in the present calcu- 
lation. For the case of Csl molecule, which we plan to take as a test molecule in the first experimental demonstration, 
the electric field of a single pulse needs to be ^ 90 kV/cm, corresponding ~ 4.3 x 10^ W/crn^ in the peak intensity. 
Under the supposed electric-field, the interaction potential of the two-photon transition is expected to be smaller than 
that of the one-photon transition by a factor of ~ lO^JLSj. The probability of the net two-photon transitions can be 
estimated from the product of the interaction potential and the number of the possible combination of the frequency 
components which satisfy the selection rule of AJ = ±2. In this case, the number of possible combination does not 
exceed J^axi— 10^). The net probability of the two-photon transitions is still smaller by a factor of ~ 10^ than that 
of the one-photon transition. We conclude that inclusion of the two-photon transition does not substantially change 
the probability distribution. 
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FIG. 4: (Color online) Simulation of the evolution induced by the optical frequency comb in the ideal situation. Exact solution 
in t = 50 is completely overlaps to the calculated lines. 



IV. COMPENSATION FOR THE CENTRIFUGAL DISTORTION 



For the case of real molecules, the pure rotational transition frequencies do not exactly compose the evenly-spaced 
comb because of the centrifugal distortion. To compensate for the difference, we have to modify the electric field defined 
in eqlH The rotational energy of the real diatomic molecules is expressed as E'j — hBJ{J + 1) — hDJ^{J + 1)^, 
where D is the centrifugal distortion constant. The transition frequency for J — J + 1 is expressed as v'j = 
2B{J + 1) - iD{J + 1)3. (For the case of Csl molecule, the ratio D/B = 1.57 x 10"'^[14].) The equation of rotation 
is re-written by 

-fije{t)cxp{-2Triiy'jt)cj+i{t). (7) 
In a similar way, the electric field for the CTQW is defined as 

e{t)^ J2 ^cos{2nv'j,t). (8) 

J'=0 

The e(t) is regarded as a train of chirped pulses whose chirp rate gradually changes from negative to positive (Fig. 5). 
Eql7]and eql8]were numerically solved by the same way as in Sec lIIII with the initial condition of cj(— 12. 5/2_B) ~ (5j.ioOi 
and D/B = 1.57 x 10~^. The calculated population distribution and the exact solution were most agreed with each 
other (Fig. 6). We conclude that the modified optical frequency comb still approximately implements the CTQW on 
the rotational states in a diatomic molecule including the effect of the centrifugal distortion. Controlling the chirp 
rate in the visible region is not difficult, and that in the mid-infrared region has been performed p^. That will be 
also possible in the terahertz region in the near future. 



V. CONCLUSIONS 



We introduced a proposal for the physical implementation of the continuous-time quantum walk using diatomic 
molecules and the optical frequency comb. We demonstrate the analogy between the CTQW and the cascade rotational 
transitions in the diatomic molecules. Possible imperfection by the centrifugal distortion could be almost eliminated 
by chirping each pulse appropriately. The implementation of the CTQW in a molecule will be useful for simulating 
the diffusion phenomena in quantum system. 
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FIG. 5: Upper panel: Time profile of the modified optical frequency comb {Jmax ~ 200, D/B = 1.57 x 10 ^, 25 pulses.) Lower 
left panel: An expanded time profile of the first pulse. Lower right panel: An expanded time profile of the last pulse. 
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FIG. 6: (Color online) Simulation of the population distribution irradiated by the modified optical frequency comb including 
the compensation of the centrifugal distortion {D/B = 1.57 x 10~^). 
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